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We develop the "simulated extinction method" to measure average foreground Galac- 
tic extinction from field galaxy number-counts and colors. The method comprises sim- 
ulating extinction in suitable reference fields by changing the isophotal detection limit. 
This procedure takes into account selection effects, in particular, the change in isophotal 
detection limit (and hence in isophotal magnitude completeness limit) with extinction, 
and the galaxy color-magnitude relation. 

We present a first application of the method to the HST WFPC2 images of the 
gamma-ray burster GRB 970228. Four different WFPC2 high-latitude fields, including 
the HDF, are used as reference to measure the average extinction towards the GRB 
in the F606W passband. From the counts, we derive an average extinction of Ay 
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= 0.5 mag, but the dispersion of 0.4 mag between the estimates from the different 
reference fields is significantly larger than can be accounted by Poisson plus clustering 
uncertainties. Although the counts differ, the average colors of the field galaxies agree 
well. The extinction implied by the average color difference between the GRB field and 
the reference galaxies is Ay = 0.6 mag, with a dispersion in the estimated extinction 
from the four reference fields of only 0.1 mag. All our estimates are in good agreement 
with the value of 0.81±0.27 mag obtained by Burstein & Heiles, and with the extinction 
of 0.78±0.12 measured by Schlegel et al. from maps of dust IR emission. However, 
the discrepancy between the widely varying counts and the very stable colors in these 
high-latitude fields is worth investigating. 

Subject headings: galaxies:photometry - galaxies:statistics - dust, extinction - ISM: 
structure 



1 Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. 
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1. INTRODUCTION 
1.1. Historical background 

Since the classical work of Trumplcr (1930a, b) pro- 
duced strong evidence of the existence of interstellar 
dust, several methods have been devised to try to 
measure the Galactic extinction towards extragalac- 
tic objects in a reliable way. In particular, there 
is a long history on the use of background galaxies 
to probe extinction. As early as 1934, Hubble real- 
ized that the variation of galaxy number-counts with 
Galactic latitude was well represented by the cosecant 
law, and was therefore consistent with a layer of ab- 
sorbing material of equal thickness above and below 
the Galactic plane. In 1958, Holmbcrg found that 
the observed colors (corrected for redshift and inter- 
nal extinction) of spirals of different types correlated 
with Galactic latitude. In 1967, Shane and Wirtanen 
published the best sample up to that date of galaxy 
number-counts; their results showed conclusively that 
galaxies were clustered, but also that Galactic extinc- 
tion was patchy. 

When compared to the colors of globular clus- 
ters and early-type field galaxies, the Shane- Wirtanen 
number-counts also raised the possibility of neutral or 
grey extinction (i.e., Rv = Ay /E(B — V) » 3), es- 
pecially close to the Galactic plane. Given that this 
"grey" extinction could be an artifact of comparing 
total extinction derived from what we cannot see with 
selective extinction (reddening) inferred from what we 
can see, Peterson (1970) used the observed colors and 
magnitudes of brightest cluster ellipticals in an at- 
tempt to measure both extinction and reddening from 
the same objects. Sandage (1973,1975) employed the 
same technique on brightest cluster members and ra- 
dio galaxies projected on the Galactic plane; he mea- 
sured Rv ~ 3, but of course the simple fact that 
he could observe the cluster galaxies indicated that 
these were regions of low extinction in the plane of 
the Galaxy. 

In a series of papers, Heiles (1976), and Burstein 
& Heiles (1978a, 1978b) discussed the relationships be- 
tween HI column density, the Shane- Wirtanen galaxy 
counts, and reddening. Heiles (1976) also addressed 
the influence of surface brightness profile and distance 
(given the dependence of apparent angular diameter 
on distance) on the probability of galaxy detection. 
Burstein & Heiles (1978b) elaborated on the effects of 
variable background when comparing number-counts; 
they also concluded that the (Poisson plus cluster- 



ing) variance in the Shane- Wirtanen number-counts 
was enough to prevent an accurate determinantion of 
extinction towards the Galactic poles. The end re- 
sult of their work was Burstein & Heiles (1982), the 
primary reference for estimating Galactic extinction 
towards extragalactic objects until the recent work of 
Schlegel, Finkbeiner & Davis (1998) using the IRAS 
and COBE/DIRBE (Diffuse Infra Red Background 
Experiment) 100 ^m maps. 

Also very recently, Gonzalez et al. (1998, hereafter 
G98) have used galaxy counts and colors to measure 
the total opacity of spiral galaxy disks. This is techni- 
cally a much more complicated problem, owing to the 
crowding and confusion from stars, star clusters, and 
HII regions in the foreground disk. In order to decou- 
ple the effects of confusion from those of extinction, 
we developed the "synthetic field" method: suitable 
and properly scaled reference fields (the Hubble Deep 
Field, HDF, in the case of G98) are added directly into 
the spiral galaxy images; the reference frames are at- 
tenuated to mimic different amounts of extinction and 
reddening, until one recovers the same number of sim- 
ulated galaxies as there are real background galaxies. 
Relative to G98, this paper is a simpler application of 
extinction simulations using reference fields. Howev- 
er, the present work also represents an advance in the 
treatment of the color information, and also because 
we use additional reference fields. 

In general, the ability to determine total average 
extinction from galaxy number-counts is limited by 
clustering, and by the difficulties in assessing com- 
pleteness. Galaxy colors, on the other hand, are bi- 
ased towards regions of low opacity (since observing 
the galaxies is a pre-condition to measure their col- 
ors). However, if it can be safely assumed that the 
region of interest can be sampled everywhere with 
field galaxies, it is possible to determine an extinction 
law from the combination of colors and counts (G98) 
or, assuming an extinction law a priori, it is feasible 
to derive total extinction from the reddening of the 
galaxies alone. This is what Zaritsky (1994) attempt- 
ed to do in the halos of external galaxies, through the 
comparison of the colors of all background galaxies 
detected with the colors of galaxies in control fields. 
But when field galaxies are thrown in together regard- 
less of type, a new problem emerges: galaxies have a 
color-magnitude dependence (Tyson 1988, Williams 
et al. 1996). This means that the color offset has to be 
measured with respect to samples that have the same 
magnitude limit as the extinction- corrected limit of 
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the galaxies detected in the field of interest. Zaritsky 
(1994) used a maximum-likelihood method to circum- 
vent this problem. We have developed a much simpler 
and more intuitive method to account for the color- 
magnitude relation when deriving Galactic extinction 
from the colors of field galaxies. 

1.2. The simulated extinction method. 

In this paper, we describe the "simulated extinc- 
tion method" to measure foreground Galactic extinc- 
tion from background field galaxy number-counts and 
colors. The method takes into account selection ef- 
fects, in particular, the change in isophotal detec- 
tion limit (and hence in isophotal magnitude com- 
pleteness limit) with extinction, and the galaxy color- 
magnitude relation. We simulate extinction in suit- 
able reference fields by changing the isophotal detec- 
tion limit. The procedure is repeated for different 
values of extinction. To make use of the cumulative 
number-counts, a plot is drawn of the total number of 
galaxies recovered in each simulation vs. the assumed 
extinction in order to ascertain the "best" fit to the 
actual number-counts in the field of interest (G98). 
To use the colors, we measure the difference between 
the average color of the galaxies in the field of inter- 
est, and the average color of the reference galaxies as 
we change the isophotal detection limit of the latter. 
Note that we are not simulating the reddening of the 
reference fields, just the change of average color with 
isophotal detection limit (the galaxy color-magnitude 
relation). Next, we plot the color difference vs. "ex- 
tinction", or detection threshold. The extinction in 
the field of interest is given by the intersection of 
this function (which decreases with extinction or a 
brighter isophotal detection limit) with the Galactic 
reddening law. 

We present a first application of the method to the 
Hubble Space Telescope (HST) Wide-Field Planetary 
Camera 2 (WFPC2) images of the gamma-ray burster 
GRB 970228. The extinction towards the GRB has 
recently been a matter of contention (Castander & 
Lamb 1998); in the absence of a measured redshift of 
the host, the assumed unreddened colors determine 
the plausibility of the object being in the Milky Way 
or extragalactic. As we will show, both galaxy counts 
and (assuming a Galactic reddening law; Schlcgel ct 
al. 1998) galaxy colors yield independently the same, 
relatively low (Ay ~ 0.6 mag) foreground extinction 
towards the burster. 



2. THE REFERENCE FIELDS 

As discussed elsewhere (G98), HST WFPC2 im- 
ages offer a significant advantage over images from 
ground-based telescopes in their greatly improved res- 
olution, which facilitates the separation of stellar ob- 
jects from galaxies. Ideally, reference fields should be 
available with the same camera and range of colors, 
and at least as deep or deeper in exposure time as the 
images of interest. Obviously, the extinction towards 
the reference fields should be known (or assumed). 
Fortunately, there are by now many fields observed 
with the WFPC2 which can satisfy these criteria, in- 
cluding the Medium Deep Survey fields (typically ex- 
posures of > 5000 s in the F606W and F8UW pass- 
bands; Griffiths et al. 1994; Ratnatunga, Griffiths, 
& Ostrander 1998), the field around the weak radio- 
galaxy 53W002 (24 orbits in 3 colors; e.g., Pascarelle 
et al. 1996a,b), the two deep Westphal fields (~ 25000 
s each in both F606W and F8UW; Westphal, Kris- 
tian, & Groth 1994) and, of course, the HDF (150 
orbits in 4 colors, Williams et al. 1996), which offers 
the deepest set of images. 

The counts and colors of the galaxies in the ref- 
erence images should be "representative" of galaxies 
in the field or at least well characterized. However, 
this is not as crucial as it might seem; as discussed 
in the conclusions to G98, the relevant quantities for 
this kind of work are the dispersions of the galaxy 
number-counts and colors, and not the means. Even 
if we manage to determine very accurately the means 
of these properties in a typical WFPC2 reference field, 
the larger source of error will still be the "single real- 
ization" of the field that we are investigating. Under 
these circumstances, the dispersion is a fair estimator 
of the expected error in the derived extinction. 

For this work, we use as reference fields the HDF, 
the field around the weak radio-galaxy 53W002, and 
the two deep Westphal fields. Like the data on the 
GRB 970228, each reference image comprises 3 con- 
tiguous Wide Field Camera fields, plus one 
Planetary Camera (PC) field which we do 
not use here Table 1 lists the equatorial co- 
ordinates, Galactic coordinates, and exposure times 
in the F606W and F8UW passbands of the refer- 
ence fields, and of the GRB 970228 images as well. 
The Burstein & Heiles (1982) E(B - V) towards all 
of the reference fields is basically zero, but we assume 

2 Actually, the optical counterpart of the gamma-ray burster was 
imaged on the PC field. 
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E(B — V) = 0.02, from the systematic offset found 
by Schlegel et al. (1998). We choose to estimate the 
clustering error from the literature (Roche et al. 1993; 
Brainerd et al. 1995), because we would need a much 
bigger sample of reference WFPC2 fields to measure 
the error directly from the data. 

3. OBSERVATIONS AND DATA REDUC- 
TION 

We use the Version 2 Release images of the HDF 
(Williams et al. 1996) retrieved from the STScI/HST 
archive. This release of the HDF was made from 
dithered exposures using the "drizzling" method 
(Fruchter & Hook 1997). The HDF was observed in 
the 4 photometric bands designated F300W, F450W, 
F606W, and F8UW (cf. the HST Data Handbook); 
we only use the images in the F606W and F814W 
passbands for this investigation (Table 1). 

Like the public release HDF images, the WFPC2 
images of the fields surrounding GRB 970228 and 
53W002, as well as the two Westphal fields, were driz- 
zled onto sub-sampled output images; care was tak- 
en to align the F606W and F814W images of each 
field. A final output pixel scale of 0705 was used, 
rather than the 0'.'04 of the HDF. (As the HDF was 
the first application of drizzling, a very conservative 
pixel scale was used.) Additionally, the images were 
cleaned of blemishes using the "blot" task (Fruchter 
& Hook 1997). Since "blot" is only used to produce 
comparison images for the creation of blemish masks, 
its use should have no effect on the photometry of the 
final images. Similarly, we are able to adjust for the 
effect of the differing pixel scales of the HDF and the 
other images by altering the source extraction soft- 
ware parameters (§4.1). 



4. ANALYSIS AND RESULTS 

When using an isophotal brightness threshold to 
detect and perform photometry of galaxies, extinc- 
tion has a double effect. It will reduce the sur- 
face brightness and the total apparent luminosity of 
the galaxies, of course, but it will also shift a fixed 
isophotal detection threshold to intrinsically brighter 
magnitudes. On average, the difference between the 
measured isophotal magnitudes and the total appar- 
ent magnitudes of galaxies will be larger than in the 
absence of extinction, because the isophotal bright- 
ness will be integrated over a smaller area of the 
galaxy (Shane & Wirtanen 1967; Heiles 1976). In oth- 



er words, the differential number-count vs. apparent 
magnitude relation will shift towards fainter magni- 
tudes by an amount larger than the extinction. An- 
other problem with comparing differential counts vs. 
magnitude in two fields is that they will be influenced 
(mainly at brighter magnitudes, where galaxy cluster- 
ing is worst) by overdensities associated with partic- 
ular large scale structures. A much better approach 
is to use cumulative counts, taking into account se- 
lection effects. 

4.1. Object extraction 

We used SExtractor (Bertin & Arnouts 1996) to 
perform the photometry on all objects in the GRB 
field, up to the 25 mag AB arcsec~ 2 isophote 
in Vf606W- This threshold is somewhat conservative, 
but it allowed the detection of most objects that we 
could identify as galaxies by eye, and at the same time 
avoided the identification of noise as spurious galax- 
ies. The source detection^, object deblending, and ob- 
ject classification (stars vs. galaxies) were done in the 
F606W image because it is deeper than the F814IF 
frame (Table 1). Also, fainter field galaxies are sta- 
tistically bluer, so searching in F606W allows the de- 
tection of more objects. Once galaxies had been iden- 
tified in the F606W data, their {F606W - F8UW) 
color was measured in a fixed aperture of 0"4 diam- 
eter (§4.3). It is worth emphasizing that separation 



between stars and galaxies is not a serious problem at 
the resolution and depth of the HST images. Not on- 
ly the separation on the basis of morphology is easier; 
in addition, owing to the finite size of the Milky Way, 
at fainter magnitudes the number-counts of Galactic 
stars are drastically reduced. We have estimated that 
the error in number-counts introduced by uncertain- 
ties in object separation is, at most, ~ 2%, about half 
of the Poisson error. 

4.2. Extinction from galaxy counts 

Following closely on the ideas developed in G98, 
photometry of galaxies was performed in the con- 
trol fields at different simulated extinctions, i.e., at 
progressively brighter detection thresholds, starting 
with the 25 Vf606W ma gAS arcsec~ 2 isophote ade- 
quate for the GRB image, and at intervals of 0.40 
mag in F606W (or 0.45 mag in V; Schlegel et al. 
1998). This will take into account the dependence 



3 A gaussian smoothing kernel with a FWHM of 2"was applied 
in all cases. 
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of the cumulative number-count vs. extinction curve 
on object selection biases. With our procedure, the 
higher signal-to-noise ratio of the reference fields (es- 
pecially the Hubble Deep Field) will result in a small- 
er photometric error, but will have no impact on the 
number of galaxies detected. We tested this assump- 
tion by adding noise to the HDF, until it had the same 
signal-to-noise ratio as the GRB image; the search on 
this degraded image yielded 277 instead of 275 galax- 
ies (Table 2). 

Figure [l] shows the comparison between the num- 
ber of galaxies found in the GRB 970228 field (hori- 
zontal line) and the number of "control" galaxies re- 
covered in the reference fields at each different sim- 
ulated extinction. The actual measurement of the 
extinction toward the gamma-ray burster field was 
performed by linearly interpolating between the sim- 
ulations with numbers of recovered galaxies immedi- 
ately above and below the number of galaxies in the 
GRB field. Panel 1 compares the galaxies in the field 
of the burster to the average of all 4 control fields; 
panels 2 through 5 show the results for the different 
reference fields. 

The error bars include the Poisson counting error 
and the galaxy clustering error expected from the two- 
point correlation function w(9 — 139")^ (Roche et al. 
1993; Brainerd et al. 1995). The correlation function 
depends on the completeness limit of the galaxies; in 
the case of the above mentioned references, w(9) is 
measured in the i?-band. We have used the HDF 
images to determine the number-count completeness 
limit to our surface brighteness detection threshold of 
25 mag^B arcsec~ 2 in F606W. From histograms of 
differential counts vs. aperture magnitude in increas- 
ingly larger apertures, we find that, in the absence 
of extinction, the number-counts are complete to a 
total integrated AB magnitude of 25.5 in F606W. 
This is equivalent to an R (apparent) brightness of 
~ 25 mag^. As we change the detection threshold 

4 The area of the three Wide Field Camera chips is equal to a 
square of 139" on the side. Since the field of view is not square, 
this slightly overestimates the value of the correlation function. 

5 We have used the IRAF/STSDAS (Tody 1993) program syn- 
phot to determine that for a K2 giant, which has (F606W — 
F81AW)ab = 0.5, that is, the same color as the average of 
all the reference galaxies to our detection limit, (F606Wab ~~ 
Rcousins) = 0.4 mag. Reassuringly, according to Smail et 
al. (1995), galaxies with a total R apparent magnitude be- 
tween 23.0 and 25.5 have (V - R) ~ 0.6 mag, or (F6Q6W A b - 
Rcouains) = 0.4. We roundoff these estimates to 0.5 mag, given 
that the references used for the two-point correlation function 



to simulate extinction in the control fields, we move 
the completeness limit accordingly, towards brighter 
magnitudes. 

Table 2 shows, for each reference field and for dif- 
ferent simulated extinctions in V, the assumed R 
magnitude completeness limit, w{9 = 139"), and the 
number of detected galaxies with their total la error; 
in parenthesis, we list separately the Poisson count- 
ing error and the galaxy clustering uncertainty. Due 
to rotation between the F606W and the F814W ob- 
servations of the GRB 970228, all the reference fields 
have an area a few percent larger than that of the 
GRB mosaic; the number-counts in Figure I and Ta- 
ble 2 have been adjusted for these differences. The 
bottom part of the table displays the number-counts 
measured in the GRB 970228 field, as well as the ex- 
tinctions in V derived from their comparison against 
the counts in the control fields. The quoted extinc- 
tions include the additional 0.07 mag in Ay implied 
by the systematic offset in E(B — V), relative to 
Burstein & Heiles (1982), found by Schlegel et al. 
(1998). Here also, we give the total error in the de- 
rived extinction and, in parenthesis, the error from 
shot-noise in the number-counts and the contribution 
to the error from galaxy clustering. 

From the average of the number-counts in all the 
control fields, we measure an extinction of Ay = 
0.51 ± 0.22 mag. We notice, however, that the dis- 
persion of the number-counts in the few WFPC2 ref- 
erence fields we use is about twice the estimated er- 
ror bars; for reasons that are yet to be understood, 
the two Westphal fields have fewer galaxies than the 
HDF and the 53W002 field, by an amount that is sig- 
nificantly larger than the expected clustering error. 
Comparison of the counts in the GRB field with those 
of the HDF and the 53W002 field gives a best esti- 
mate of the extinction of Ay — 0.75±0.24; if the two 
Westphal fields are averaged in with the richer fields, 
and we use the dispersion to estimate the error, the 
best value of the extinction is Ay — 0.51 ± 0.37 mag. 

4.3. Extinction from galaxy colors 

We measured the average colors of the reference 
galaxies, once again at progressively brighter detec- 
tion thresholds, starting with the 25 Vp606VK ma gAB 
arcsec -2 isophote adequate for the GRB image. Fol- 
lowing G98, after identification in the F606W images, 



provide their measurements in half-magnitude steps. If any- 
thing, this slightly overestimates the expected clustering error. 
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we measured the color of each individual galaxy in a 
fixed aperture, in this case of 0"4 diameter. The goal 
is to avoid contaminating the measured galaxy col- 
ors with the color of the sky. This means that the 
main effect of the change in detection threshold (and 
therefore in isophotal magnitude limit) is the loss of 
faint (statistically bluer) galaxies. The error intro- 
duced by not measuring possible color gradients over 
larger galaxy radii should be negligible, especially be- 
cause the measurement of a color gradient depends 
very critically on an accurate estimation of the back- 
ground 0. We assume that the error in the mean color 
is its measured dispersion divided by y/N, where N 
is the number of galaxies used for the color measure- 
ment []. 

Since the errors in the average colors are of the 
order of 0.05 mag (Table 3), in the simulations we 
incremented the detection threshold in steps of 0.10 
mag in F606W, or 0.11 mag in V (Schlegel et al. 
1998). We remind the reader that we did not simulate 
the reddening of the reference fields, just the change 
of average color with isophotal detection limit (the 
galaxy color-magnitude relation) . We then calculated 
the color difference (A Color) in (^606^-^8141^) 
between the reference galaxies and the galaxies in the 
GRB field. 

Figure^ shows the change in A Color with isopho- 
tal detection limit. Panel 1 compares the galaxies in 
the field of the burster to the average of all 4 control 
fields; panels 2 through 5 show the results for each 
one of the different reference fields. A Color decreas- 
es at brighter detection thresholds (i.e., higher extinc- 
tion), due to the color-magnitude relation shown by 
field galaxies (brighter galaxies are statistically red- 
der). The errors are higher at high extinction, both 
because fewer galaxies are detected and because the 
intrinsic dispersion in color is higher for apparently 
brighter galaxies. Figure || also displays the Galac- 
tic reddening line, that is, E(F606W - F8UW) vs. 
Ay (Schlegel et al. 1998). The intersection between 
the two relations measures the foreground Galactic 

6 When using isophotal apertures, the average color of the galax- 
ies in each field changes by less than the photometric error, 
relative to the average color measured in the 0"4 apertures. 

7 Given that the color distribution of galaxies is not gaussian, 
in G98 we compared the gaussian estimate of the error to the 
dispersion in the mean color of 100 random samples of N HDF 
galaxies (where N is the number of galaxies revcovered at each 
simulated extinction); we found that the error predicted from 
bootstrapping statistics is of the same order as the error found 
by assuming a gaussian distribution. 



extinction A v towards GRB 970228. The observed 
average color of the galaxies in the GRB field equals 
the sum of the control field average color plus the red- 
dening at the inferred extinction (cf. Table 3). Albeit 
much simpler, our procedure is equivalent to simul- 
taneously changing the isophotal detection threshold 
and reddening the control galaxies following a Galac- 
tic extinction law, until their average color is the same 
as that of the galaxies in the GRB field. This is so be- 
cause we base the galaxy detection in only one band, 
and we do not use an upper limit in color as a galaxy 
selection criterion. 

Table 3 shows, for each reference field and for dif- 
ferent simulated extinctions in V, the average color of 
the galaxies and the difference in color with respect 
to the galaxies in the GRB field. The bottom of the 
table displays the measured color of the galaxies in 
the GRB field, and the derived extinctions in V from 
the comparison with the galaxy colors in the refer- 
ence fields. Accounting as described for the color- 
magnitude relation of field galaxies, and assuming a 
foreground screen Galactic reddening law, a compar- 
ison between the average colors of, respectively, the 
control galaxies and the galaxies in the GRB field 
yields Ay — 0.61±0.11 mag; the error represents the 
scatter between the measurements obtained from the 
different control fields. In contrast to the number- 
counts, the galaxy colors in all 4 reference fields agree 
well; within the errors, the colors of all the fields are 
the same, at all detection thresholds. 

5. SUMMARY 

We have developed the "simulated extinction 
method" to measure foreground Galactic extinction 
from the number-counts and colors of field galaxies. 
The method simulates extinction in suitable reference 
fields by changing accordingly the isophotal detection 
limit. This procedure allows the calibration of selec- 
tion effects and systematics, in particular the galaxy 
color-magnitude relation, and the change in isophotal 
detection limit with extinction (and hence the change 
in both the isophotal magnitudes of galaxies and in 
the number-count completeness limit). 

We have used 4 HST WFPC2 reference fields to 
measure forward Galactic extinction towards the field 
of the gamma-ray burster GRB 970228. Compari- 
son of the counts in the GRB field with those of the 
HDF and the 53W002 field gives a best estimate of 
the extinction of A v = 0.75±0.24; if the two West- 
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phal fields arc included, the value of the extinction is 
Ay — 0.51 ± 0.37 mag. Although the counts differ, 
the average colors of the field galaxies agree well. The 
extinction implied by the average color difference be- 
tween the GRB field and the 4 reference fields is Ay 
= 0.61 mag; the dispersion in estimated extinction 
between the 4 reference fields is only 0.11 mag. How- 
ever, there may be a slight bias in the value measured 
from the colors; if the extinction is clumpy, the most 
obscured galaxies may fall out of the sample and not 
affect the average color. Thus the estimate of Ay = 
0.61 mag may be a lower limit. 

6. CONCLUSIONS AND DISCUSSION 

All our measurements are in good agreement with 
the value of 0.81±0.27 mag obtained by Burstein & 
Heiles (1982; error re-estimated by Schlegel et al. 
1998), and with the extinction of 0.78±0.12 mag mea- 
sured from maps of dust IR emission (Schlegel et al. 
1998). Ultimately, all of these estimates are averages 
with different spatial resolutions and over different 
components of the ISM, at a low latitude position 
where the Milky Way exhibits a steep gradient in ex- 
tinction (cf. Burstein & Heiles 1982); none of these 
methods is capable of measuring the extinction in the 
exact line of sight to the GRB. 

It is very exciting that the colors of the field galax- 
ies agree so well, despite the big differences in number- 
counts and even though, for example, there are no 
clusters at low redshift in the HDF, while the 53W002 
field was selected around a weak radio galaxy (in no 
case did we try to identify and excise galaxy clusters). 
Moreover, comparison of the extinction-corrected col- 
or of the GRB host to the {F606W - H) colors of the 
HDF galaxies (Fruchter et al. 1998) suggests that the 
value of Ay = 0.6 mag derived from the colors might 
be closest to the correct extinction; already with the 
value of Ay = 0.75 mag adopted in that paper, the 
GRB host is as blue as the bluest objects in the HDF. 

But although the tightness in the colors is very in- 
teresting per se and from the point of view of their 
usefulness to probe extinction, the discrepancy be- 
tween the widely varying counts and the very stable 
colors in these high-latitude fields is worth investigat- 
ing. 

The 53W002 field has a slightly lower exposure 
time than the two Westphal fields. This, and the 
galaxy search performed in the degraded version of 
the HDF (§4.2) eliminate the possibility that the dis- 



crepancy is caused by the different signal-to-noise ra- 
tios of the reference fields. The disparity occurs across 
all magnitude bins fainter than 21.5 Vfgo6W niag 
(isophotal), and applying a fainter detection thresh- 
old does not really help to reduce it. We looked at the 
counts and colors of the reference galaxies when us- 
ing a surface brightness detection limit of 26 mag^ 
arcsec -2 isophote in Vf&oqw (a threshold set by the 
53W002 field, i.e., the shallowest of the reference 
images). The colors remain as stable; the number- 
counts, on the other hand, start to converge, but 
not nearly as quickly as the expected clustering er- 
ror drops. 

It is possible that the problem lies in the quot- 
ed clustering errors. To our knowledge, the field-to- 
field variance on the scale of the WFPC2 field of view 
has not been investigated, even on ground-based da- 
ta. In fact, the characteristics of fields of the size, 
spatial resolution, and depth of our reference images 
are a matter of intense study at present. In this re- 
gard, although the properties of the HDF (Villumsen, 
Freudling, & Da Costa 1997) appear to be consistent 
with those derived for ground-based data (Brainerd 
et al. 1995), results on other deep fields are not yet 
available. The clustering of galaxies brighter than I 
= 23 mag in prerefurbishment WFPC Medium Deep 
Survey images (Neuschaefer et al. 1995) is also consis- 
tent with ground-based data (Roche et al. 1993), but 
the completeness limit of the survey is much brighter 
than that of the WFPC 2 deep fields. 

We could also have been too optimistic in our es- 
timate of the completeness limit. If we disregard the 
results of the aperture photometry and instead take 
the view of Tyson (1988), that the completeness lim- 
it of a survey is about two magnitudes brighter than 
the surface brightness detection threshold, then our 
number-counts are complete to only R ~ 22.5 in the 
absence of extinction. This will in fact bring our er- 
rors in line with the observed dispersion of the data. 

Lastly, 4 fields are not that many. We cannot rule 
out the possibility that we have been extremely un- 
lucky and are just observing a large statistical fluctu- 
ation. 

It is important to find out whether the Westphal 
fields are outlyers or, conversely, if the observed dis- 
persion in the number-counts reflects the real disper- 
sion between WFPC2 fields. If the latter is true, the 
size of the errors severely hampers the precision of ex- 
tinction measurements derived from galaxy number- 
counts. Also, studies of the clustering properties of 
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a number of deep WFPC2 fields should determine 
whether the inadequacy of the estimated clustering 
errors is real or caused by an overly optimistic com- 
pleteness limit. Finally, even if it sounds improbable, 
there is of course the possibility that some obscuring 
material lies in between us and the Westphal fields. It 
would have to be colder than ~ 10-15 K to not show 
up in the 100 fim maps (Schlcgcl et al. 1998), but this 
is not unlikely at high Galactic latitude, far from the 
interstellar UV radiation field; it would also have to 
leave the colors of the galaxies unchanged. The an- 
swer to the number-count discrepancy is beyond the 
scope of this paper, but it will be very interesting, 
either way. 

We acknowledge Max Mutchler for assistance with 
the reduction of the GRB WFPC2 images. We are 
indebted to Ron Allen and Harry Ferguson for their 
very useful comments. We thank Dave Schlegel for 
providing measurements of the reddening towards all 
the fields we use in this paper, and Paul Goudfrooij for 
help with synphot. Support for this work was provid- 
ed partly by NASA through grant AR-06400.01-95A 
and partly by the Director's Discretionary Research 
Fund at the Space Telescope Science Institute. 
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Table 1 
Observation Log. 















Exposure time (s) 


Field 


a (J2000) 


6 (J2000) 


I 


b 


F606W 


F8UW 


GRB 970228 


05h 01m 46.7s 


11° 46' 


53" 


188.91° 


-17.94° 


9400 


4800 


HDF 


12h 36m 49.4s 


62° 12' 


58" 


125.89° 


54.83° 


109050 


123600 


53W002 


17h 14m 14.9s 


50° 15' 


30" 


77.02° 


35.79° 


20400 


20400 


Westphal 1 


15h 58m 49.8s 


42° 05' 


23" 


66.93° 


49.23° 


25200 


25200 


Westphal 2 


14h 17m 43.6s 


52° 28' 


41" 


96.35° 


59.97° 


24400 


25200 
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Table 2 
Galaxy Counts 



Reference Fields. 



Input A v * 


7?. limit 


w{6 = 139") 






N 












< All > 


HDF 


53W002 


Westphal 1 


Westphal 2 


0.07 


~ 25.0 


0.006 b 


202 ± 14 (7,12) 


275 ± 36 (16,32) 


229 ± 30 (15,26) 


150 ± 21 (12,17) 


156 ± 21 (13,17) 


0.52 


~ 24.5 


0.007 b 


143 ± 11 (6, 9) 


200 ± 29 (14,25) 


164 ± 24 (13,20) 


97 ± 15 (10,11) 


113 ± 17 (11,13) 


0.97 


~ 24.0 


0.010 d 


89 ± 8 (5, 7) 


124 ± 22 (11,19) 


104 ± 18 (10,15) 


55 ± 11 ( 7, 8) 


71 ± 13 ( 8,10) 


1.42 


~ 23.5 


0.012 c 


58 ± 6 (4, 5) 


74 ± 15 ( 9,12) 


72 ± 14 ( 8,12) 


34 ± 8 ( 6, 6) 


53 ± 11 ( 7, 9) 


1.87 


~ 23.0 


0.020 c 


35 ± 5 (3, 4) 


47 ± 12 ( 7,10) 


42 ± 11 ( 6, 9) 


18 ± 5 ( 4, 4) 


32 ± 9 ( 6, 7) 


GRB 970228. 




7?. limit 


w(0 = 139") 






Derived A v * 












All 


HDF 


53W002 


Westphal 1 


Westphal 2 


145 ± 22 (12,18) 


~ 24.5 


0.007 b 


0.51 ± 0.22 f (0.12,0.19) 


0.84 ± 0.34 (0.17,0.29) 


0.66 ± 0.34 (0.18,0.29) 


< 0.42 e (0.18,0.25) 


< 0.49 e (0.17,0.25) 



COLS.— Top: Input attenuation (mag) in the V-band; (2) 7?-band completeness limit (apparent mag); (3) two-point correlation function; (4) - (8) number-counts 
with total la error and, in parenthesis, counting Poisson error and galaxy clustering uncertainty. Bottom: (1) number-counts; (4) - (8) derived attenuation with total 
lcr error and, in parenthesis, error from the shot-noise in the number-counts and error contributed by galaxy clustering. 

a All values of A v include offset of +0.02 mag in E(B - V) from Schlegel et al. (1998). 

b Brainerd et al. (1995). 

c Roche et al. (1993). 

interpolated between values in the two previous references. 
e Value already includes +la = 0.31 mag. 

f Error includes Poisson uncertainty and theoretical clustering; the dispersion of the extinctions derived from the 4 reference fields is 0.37 mag. 



Table 3 

Galaxy <F606W - F814W> Colors. 



Reference Fields. 



Input Av a Reddening A Color Color 











All 






HDF 


53W002 


Westphal 1 


Westphal 2 


All 


HDF 


53W002 


Westphal 1 


Westphal 2 


0.07 


0. 


00 


0. 


22±0 


04 


0. 


23±0.05 


0.23±0.05 


0.23±0.05 


0.15±0.05 


0.50±0.01 


0.49±0.03 


0.49±0.03 


0.46±0.03 


0.55±0.03 


0.18 





04 


0. 


20±0 


04 


0. 


22±0.05 


0.22±0.05 


0.23±0.05 


0.14±0.05 


0.51±0.02 


0.50±0.03 


0.50±0.03 


0.48±0.03 


0.56±0.03 


0.29 





.07 





19±0. 


04 


0. 


21±0.05 


0.21±0.05 


0.22±0.05 


0.13±0.05 


0.52±0.02 


0.49±0.03 


0.50±0.03 


0.50±0.03 


0.58±0.04 


0.40 
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18±0. 


.04 


0. 


20±0.05 


0.19±0.05 


0.21±0.05 


0.12±0.06 


0.53±0.02 


0.51±0.03 


0.50±0.03 


0.51±0.04 


0.61±0.04 


0.51 





.13 





16±0. 


04 


0. 


18±0.05 


0.18±0.05 


0.20±0.06 


0.11±0.06 


0.55±0.02 


0.53±0.03 


0.52±0.03 


0.54±0.04 


0.63±0.04 


0.63 





.17 


0. 


15±0. 


04 


0. 


17±0.05 


0.16±0.05 


0.19±0.06 


0.10±0.06 


0.56±0.02 


0.53±0.03 


0.55±0.03 


0.54±0.04 


0.62±0.04 


0.74 





20 





15±0. 


.05 


0. 


16±0.05 


0.15±0.05 


0.18±0.06 


0.09±0.06 


0.56±0.02 


0.55±0.03 


0.57±0.03 


0.52±0.04 


0.62±0.05 


0.86 





.23 





14±0. 


.05 


0. 


15±0.06 


0.13±0.05 


0.17±0.06 


0.08±0.06 


0.57±0.02 


0.55±0.04 


0.58±0.04 


0.54±0.05 


0.64±0.05 


0.97 





.27 


0. 


12±0. 


05 


0. 


14±0.06 


0.12±0.06 


0.16±0.06 


0.06±0.07 


0.59±0.02 


0.59±0.04 


0.58±0.04 


0.55±0.05 


0.66±0.05 


1.08 





30 


0. 


11±0. 


05 


0. 


12±0.06 


0.10±0.06 


0.16±0.06 


0.05±0.07 


0.61±0.02 


0.60±0.04 


0.62±0.04 


0.54±0.05 


0.66±0.06 


1.20 





33 





09±0. 


05 


0. 


11±0.06 


0.09±0.06 


0.15±0.06 


0.04±0.07 


0.62±0.02 


0.61±0.04 


0.64±0.04 


0.56±0.05 


0.63±0.06 



GRB 970228. 


Color" 


All 


HDF 


Derived Ay*" 
531^002 


Westphal 1 


Westphal 2 


0.71±0.04 


0.61±0.11 d 


0.65±0.06 


0.63±0.05 


0.71±0.08 


0.46±0.05 



Cols.- Top: (1) input attenuation (mag) in the V-band; (2) assumed reddening E(F60(>W - F81AW) for Galactic extinction; (3) - (7) difference (mag) between the 
average {F&mW - F8UW) color of control galaxies and of galaxies in the GRB 970228 field: (8) - (12) average (FGOdW - F8UW) color of galaxies in the different 
reference fields. Bottom: (1) average color of galaxies in the GRB 970228 field; (2) - (6) Galactic extinction (mag) towards the GRB, derived from intersection of 
A Color with the Galactic reddening curve (see text and Figure 2). 

a All values of A v include offset of +0.02 mag in E(B - V) from Schlegel et al. (1998). 
b Schlegel et al. (1998). 

c The observed average color of the galaxies in the GRB field equals the sum of the control field average color plus the reddening at the inferred extinction. 
d Error represents scatter between the 4 reference fields. 
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Fig. 1. — Average number of recovered galaxies vs. simulated extinction of Ay- Top left: All control fields; top 
right: HDF; middle left: 53W002; middle right: Westphal 1; bottom left: Westphal 2. The error bars include the 
statistical Poisson uncertainties, as well as the field-to-field variations expected from galaxy clustering. Horizontal 
solid line: number of galaxies in the GRB 970228 field, with its uncertainty (dashed lines). For fixed number-counts, 
the error from galaxy clustering would increase at higher extinctions because the completeness limit would move 
towards intrinsically brighter magnitudes, where clustering is worse. 
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Fig. 2. — Difference in average color between reference galaxies and galaxies in the GRB field. Top left: All control 
fields; top middle: HDF; top right: 53W002; bottom left: Westphal 1; bottom middle: Westphal 2. The difference 
diminishes as the detection threshold moves to brighter magnitudes due to extinction; the error bars assume a 
gaussian distribution of the colors. Dotted line: Linear fit to the color difference; solid line: Galactic reddening law 
(Schlegel et al. 1998). The intersection of these two lines yields the extinction towards the burster (see text). 
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